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The scaling of the decoherence factor of a qubit coupled to a spin chain driven across 

quantum critical points 

Tanay Nag, 1 Uma Divakaran, 1 and Amit Dutta 1 
department of Physics, Indian Institute of Technology Kanpur, Kanpur 208 016, India 

We study the scaling of the decoherence factor of a qubit (spin— 1/2) using the central spin model 
in which the central spin (qubit) is globally coupled to a transverse XY spin chain. The aim here is 
to study the non-equilibrium generation of decoherence when the spin chain is driven across (along) 
quantum critical points (lines) and derive the scaling of the decoherence factor in terms of the driving 
rate and some of the exponents associated with the quantum critical points. Our studies show that 
the scaling of logarithm of decoherence factor is identical to that of the defect density in the final 
state of the spin chain following a quench across isolated quantum critical points for both linear 
and non-linear variations of a parameter even if the defect density may not satisfy the standard 
Kibble-Zurek scaling. However, one finds an interesting deviation when the spin chain is driven 
along a critical line. Our analytical predictions are in complete agreement with numerical results. 
Our study, though limited to integrable two-level systems, points to the existence of a universality 
in the scaling of the decoherence factor which is not necessarily identical to the scaling of the defect 
density. 

PACS numbers: 



When a quantum many-body system is slowly driven 
across a quantum critical point (QCP)i by varying a 
parameter in the Hamiltonian, defects are generated in 
the final state; this is a consequence of the diverging 
relaxation time close to the QCP, so that the dynam- 
ics is no longer adiabatic however slow may the vari- 
ation be 2 - - — . If a parameter A of the Hamiltonian de- 
scribing a d— dimensional system is changed linearly as 
\{t) = t/r, -oo < t < oo, (with the QCP at A = 0), the 
defect density (n) in the final state satisfies the Kibblc- 



Zurck (KZ) scaling relation 
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r is the inverse rate of quenching, and v and z are the 
correlation length and dynamical critical exponents, re- 
spectively, associated with the QCP. 

In parallel, there are a plethora of studies which con- 
nect quantum information theory to quantum critical 
systems (for a review, see^). One of the major issues 
in this regard is the study of decoherence, namely, the 
loss of coherence in a quantum system due to its interac- 
tion with the environment^. To elucidate these studies, 
the central spin model (CSM) has been proposed— In 
this model, a central spin (CS) (i.e., the qubit) has a 
global interaction with a quantum many body system 
(e.g., with all the spins of a quantum spin chain) which 
acts as the environment. The interaction between the 
qubit and the environment in fact provides two channels 
of time evolution of the environmental spin chain. It has 
been observed that the purity of the CS is given in terms 
of the Loschmidt echo (LE) or the decoherence factor 
(DF) which is the measure of the square of the overlap of 
the wave function evolved along the two different chan- 
nels as a function of time. The LE or DF which appears 
in the off-diagonal term of the reduced density matrix 
of the qubit is minimum at the QCP signifying a maxi- 
mum loss of coherence close to it which can be used as 
an indicator of quantum criticality-i^r— . 

At this point, the natural question would be what hap- 



pens when the environment is driven following some pro- 
tocol across a QCP. In a recent work, Damski et ati&, 
studied the decoherence of the CS by coupling it to a 
transverse Ising spin chain which is driven across the 
QCP by a linear variation of the transverse field and 
showed that in the limit of weak coupling, the logarithm 
of the non-adiabatic part of the DF (arising due to the 
contribution of the low-energy modes close to the criti- 
cal mode and denoted by D non - ac i) satisfies an identical 
scaling to that of n, given by In D non ^ a d ~ t -1 / 2 . 

In present work, we consider a version of the CSM 
in which environmental spin chain is chosen to be an 
anisotropic XY spin chain which has a rich phase di- 
agram and thereby enables us to study the scaling of 
the DF when the environment is quenched across dif- 
ferent critical and multicritical points as well as gaplcss 
critical lines. The spin chain is exactly solvable using 
Jordan- Wigncr(JW) transformations^ and is reducible 
to a decoupled two-level problem in the Fourier space 
and hence the Schrodinger equations describing the evo- 
lution of these two levels can be analytically solved for 
all times. We shall, however, emphasize on an alternative 
method introduced in— This method, valid away from 
the QCP, exploits the two-level nature of the reduced 
Hamiltonian and the exact expression for the probability 
of non-adiabatic transition at the final time as given by 
the Landau-Zener (LZ) transition formula—. Both the 
results however lead to identical scaling relations which 
are also verified numerically. Moreover, this alternative 
approach also allows us to calculate the scaling of the 
DF for a non-linear variation of the quenching parame- 
ter though the LZ formula is not known exactly. 

Let us emphasize that our focus here is limited to only 
low-energy modes close to the critical mode for which the 
energy gap vanishes at the QCP. The high energy modes 
on the other hand, evolve adiabatically throughout the 
dynamics; though these modes contribute to the dynam- 



ics of dccoherence non-trivially through the fidelity fac- 
tor, they do not alter the scaling relation of the DF— . 

Let us first clarify the connection between the scal- 
ing of the DF and n, that we are interested in. We 
shall assume weak coupling between the qubit and the 
environment and work within the appropriate range of 
time; under these circumstances, for all the quench- 
ing schemes discussed here (achieved by changing a pa- 
rameter A = t/r), we find the scaling relations: (i) 
ln£)„ „_ ad ~ {-t 2 f{r)), if QCP is at A = and (ii) 
lnL>„on-ad ~ {-(* - A r) 2 /(r)}, if the QCP is at A 
(as happens for quenching through a MCP discussed 
below). We explore the scaling of this function /(r) 
(which is found to be linear in the size of the spin 
chain and quadratic in system-environment coupling) 
with r and address the question whether that is iden- 
tical to the scaling of n. However, to eliminate t, one 
could further substitute t = At, to obtain the scaling 

\nD non _ ad A 2 t 2 /(t) (or (A- A ) 2 t 2 /(t) for case 

(ii)), but it should be emphasized that the non-trivial 
scaling of lnD no „_ a ^ with r is provided by that of /(r). 
Our studies reveal that in the cases when the environ- 
ment is driven through an isolated QCP or a multicriti- 
cal point (MCP), the scaling of \nD non ^ a d (or precisely 
that of /(t)) is the same as that of the defect density 
in the final state following a quench for both linear and 
non-linear quenches. However, there are situations when 
this generic connection do not hold. For example, when 
the environment is driven along a critical line across the 
MCP, we arrive at a scaling which is significantly differ- 
ent from that of n. 

The Hamiltonian He of the environment is the XY 
spin chain in a transverse field consisting of N spins given 
byii 
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H e = - ]T[ J -«+i + J y afaf +1 + hat] 
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which is coupled to the spin- 1/2 qubit by a Hamiltonian 
Hse- In the following, we shall define the anisotropy 
parameter 7 = J x — J y , and the parameter (J x + J y ) will 
be set equal to unity in all cases except for the quenching 
through a MCP. The model is exactly solvable by JW 
transformation^; the phase diagram is shown in FigQ] 

Let us introduce the notion of DF by considering the 
situation in which the transverse field h is quenched as 
h(t) = 1 — t/r, and the qubit is coupled to the time de- 
pendent transverse field of ((TJ) through the Hamiltonian 

Hse = — S^2i = i a i a si where °f is the i— th spin of the 
XY chain and a s represents that of the qubit. The sys- 
tem crosses the Ising critical points at h = 1 and h = — 1 
with critical mode k c given by k c = n, and 0, respectively. 
We choose the qubit to be initially (at t — > —00) in a pure 
state superposition \4>s{t — > —00)) = ci| t)+ c 2| I), where 
I f) an d I |) represent up and down states of the CS, re- 
spectively, and the environment is in the ground state 
\<pE(t —> —00)) = \<f> g ). The ground state of the compos- 
ite Hamiltonian He + Hse, at t — > —00, is given by the 
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FIG. 1: (Color Online) The phase diagram of Hamiltonian JT]) 
in the h/(J x + J y ) — j/(Jx + Jy) plane. The two red vertical 
lines are the transverse Ising critical lines between the fer- 
romagnetic (FM) and paramagnetic (PM) phases. The blue 
horizontal line for — 1 < h < 1 corresponds to the anisotropic 
critical line separating the two ferromagnetic phases FM X and 
FM y , with ordering in x and y directions, respectively. Points 
denoted by A and B represent the two multicritical points. 
Green lines represent different quenching paths discussed in 
this paper. 



direct product \ij)(t —¥ —00)) = \<f>s(t —> —00)) ® \4>g)- It 
can be shown that at a later time t, the composite wave 
function is given by \ip(t)) = a\ f) ® |</>+) +c 2 | I) <g) \<j>-), 
where \<f>±) are the wavefunctions evolving with the envi- 
ronment Hamiltonian i7g(ft,±(5) given by the Schrddinger 
equation id/dt\(j>±) = H[h(t) ± <5]|0±). We therefore find 
that the coupling 6 essentially provides two channels of 
evolution of the environmental wave function with the 
transverse field h + 6 and h — 5, respectively. 

It is straightforward to show that the decoherence fac- 
tor D{t) defined as \(4> + (t)\4>-{t)}\ 2 , is the off-diagonal 
element of the reduced density matrix of the qubit^. To 
evaluate D(t), we rewrite the Hamiltonian ((T|) with mod- 
ified h (due to the coupling 5) in terms of JW fermions 
which then can be decoupled into a sum of independent 
(2 x 2) Hamiltonians in the Fourier spac o 14 ' 17 . In the ba- 
sis |0) and \k, —k), which represent no quasiparticle, and 
quasiparticlcs with momentum k and —k, respectively, 
the Hamiltonian He can be written as 



H%{t) = ^H±(t), where, 



HHt) 



( h{t) ±5 + cosh 






7 sm / 



7sin/c 
-{h{t)± 8 + cos k) 



The general wave function for He at any instant t can 
be written as 

I^W) = IH^(*)> = II [«*C*)|o> + «* (*)!*»-*>] - 

k fe>0 

The coefficients u^ and vj~ are obtained by solv- 
ing the Schrddinger equation id/dt(u^(t),v^(t)) = 

Hk{t)(ufc(t),Vfc(i)) where A T represents the trans- 
pose operation of the row matrix A. Hence, 



the expression of D(t) is given by rifc-^feM 
U k \(MHt) + S)\Mh(t)-S))\ 2 ,ov, 



D(t) = exp 



N 

2^ 



dk lnF, 



(2) 



where F k can be written in terms of u k and v k . We re- 
iterate that we shall focus in the limit of small 5 and 
consider only the low-energy modes which show non- 
adiabatic behavior close to the QCP. On the other hand, 
the high energy modes evolve adiabatically and their 
overlap is close to unity. This method can be useful for 
exact solution as well as numerical estimation of D(t). 

We shall however introduce a simpler method for an- 
alytical calculations that exploits the (2 x 2) nature 
of the reduced Hamiltonian to calculate i^ (£)!£. Far- 
away from the QCP (|/i(£)|»l (t — > +oo)) i.e., after 
crossing both the QCPs, we can write \4>k(h + 5)) = 
Ufc|0) + v k e- iA+t \k,-k), and \< j> k {h - 5)) = u k \0) + 
e~ iA ~'w fe |fc,-fc) where A+ = 4y/(h + S + l) 2 + j 2 sin k 2 
and A - = kyih — 8 + l) 2 + 7 2 sin k 2 are the energy dif- 
ference between the states |0) and \k,—k) when the trans- 
verse field is equal to h + S and h — 5, respectively. In 
writing the above expression, we make use of the fact that 
excitations occur only in the vicinity of QCPs. Follow- 
ing that the wavefunctions (|0 ± (i))) evolve adiabatically 
picking up the appropriate phase factor with time. At 
the same time, the coefficients u k and v k can be found to 
be |iifc| 2 = 1 — Pk and |t>fc| 2 = Pk where pk is the Landau- 
Zener probability of excitations for the mode k given by 
Pk = cxp(— 2-7TT7 2 sin 2 k)^. Combining all these, we find 



F k {t) = \(MHt) + s)\MHt) - s))f 



Kr + kre 



2 p -i(A+-A-)t 



(3) 



which can be recast in the vicinity of the quantum critical 
point at ft = 1 to the form A = (A+ - A _ )/2, 

F k (t) = l-4p k {l-p k )sm 2 {At) 



= 1-4 



n 2,/2 

-2-7TT7 K 



-47TT7 k ' 



sin 2 (4<«) (4) 



where sin k has been expanded near the critical modes 
fc = 7r, with k' = n — k and we have taken the limit 
6 — > 0. The above expression is identical to that given 
ir*i^ derived via the exact solution of the Schrodinger 
equation. 

The DF is the product of the contribution from the 
modes evolving adiabatically (given by fidelity) and the 
modes evolving non-adiabatically denoted byD ncm ^ ad . 
The expression of D non - ad (t) due to the non-adiabatic 
dynamics of modes k ~ 7r after crossing the critical point 
h = 1 can be obtained from Eq. ([2]) in the following way: 
in the limit <5 — s- (or more precisely (St) — > 0), one can 
approximate sin 2 4St sa 16S 2 t 2 which results to 



N f°° 

Dnon-adtt) = exp — - / dk 

2lT Jo 



In 1 — ' P - 2lTT i" k '' — P -^ T 



^ 2fc ' 2 ) 64<5 2 i 2 }5) 



where we have extended the limit of integration to oo 
since only the modes close to the critical modes con- 
tribute in the limit of large r. Using the fact that 
ln(l — x) ~ —x, for small x, it can be further shown 
that D non _ ad is given by 



D non - ad (t) - cxp{-8(V2 - l)N5 2 t 2 /(jir^)}- 



(6) 



It is worth noting that the periodicity in time as in 
Eq. (d|) is lost and there is an exponential decay as shown 
in Eq. (j6)). This Gaussian form holds true when t <C 
1/(5 and 5 — > 0; clearly the time range over which this 
is applicable increases with decreasing 6. Otherwise, a 
sinusoidal variation is observed. A similar expression can 
be obtained for the low k modes excited after crossing the 
h = —1 critical point. We find that lnD non _ ad (t) rather 
/(t) = {8(V2- 1)NS 2 }/ (it 7 </¥) varies as l/v 7 ^, a scaling 
which is identical to that of n (with d = v = z = 1). 

Above calculations can be extended to the case of 
the non-linear quenching of a term of the Hamiltonian, 
e.g., with the variation of the transverse field h given 
by 1 — sgn(t)(t/r) a , where sgn stands for the sign func- 
tion of t. Although the probability of excitation is not 
exactly known, casting the Schrodinger equations which 
describe the time evolution of the two-level systems to a 
dimensionless formic, it has been argued that p k should 
be a function of the dimensionless combination of k and 
r given by pk = G(fc 2 r 2Q /( Q+1 '), where G is the scal- 
ing function. Considering only the contributions from 
the low-energy modes for large r, one finds the seal- 
ing D non . ad {t) = exp(-CN8 2 t 2 /T a /^+^), where C is 
a number which also depends on a. This is again in con- 
gruence with the scaling of n for a non-linear quenching 
i.e., n ~ T~ a /( a+1 \ This scaling has been numerically 
verified by directly integrating the Schrodinger equation 
(see discussion around Eq. ([2])) and results are presented 
in Fig. (Ha). 

In order to extract the exponent in a transparent way 
from the numerical data, double logarithm of D non - ac [ 
is required which is numerically not possible since D is 
always less than unity. Hence, to calculate the expo- 
nent of r, we introduce a modified DF, A(t, t) given by 



A(r,t) 



logio D 



i — ad- 



Fig. [2ja) clearly shows that 



In A(t, t) varies linearly with In r and has a slope given 
by —a/(a + 1) with a fixed t, thus confirming the ana- 
lytically predicted scaling relation. 

If the parameter h is set equal to 2 J y and the inter- 
action term J x is quenched as t/r, the spin chain (TT]) is 
driven across the quantum MCP A at J x — J y or t = J y r 



(see Fig. []}; n satisfies a scaling relation^ n ~ r 



-l/6_ 

This is not in agreement with the KZ prediction and has 
been justified by asserting the existence of quasi-critical 
points on the ferromagnetic side of the MCP—. What 
would happen if the environmental spin chain is driven 
across the MCP? Choosing appropriately the interac- 
tion Hsb^-, one finds lnD non - ad (t) ~ (t- J v t) 2 /t 1/6 ~ 
(Jx-J s )r 11/6 (see Fig. Hb)). 

Our studies have so far been limited to isolated quan- 
tum critical and multicritical points and in all cases, the 




FIG. 2: (Color Online) (a) The variation of In A(t) with lnr 
for two different values of a with N = 300 and <5 = 0.0001: 
(i) a = 0.8 (red line) for fixed t = 1500 with slope -0.42 
and, (ii) a = 1.2 (green dashed line) for t — 670 with slope 
—0.54; both the slope values are very close to the predicted 
—a/ (a + 1). (b) The same for three different types of qubit- 
environment coupling with 8 = 0.0001 and a — 1. (i) Line 
(red) corresponds to quenching across the multicritical point 
A in Fig. [l]by varying J x ~ t/r with J x — 7 and N = 500; the 
slope ~ 1.9 (ii) Dotted line (pink) corresponds to quenching 
7 ~ t/r at h — 0.5 (across the anisotropic QCP) with 7 = 6.5 
and N=200, and the slope ~ 1.4 (iii)Dashed line (green) is 
for 7 quenching as t/r along the gapless line h — 1 with slope 
~ 0.9, for 7 = 5 and N=400. Analytically predicted values 
are 11/6, 3/2 and 1, respectively. 



scaling of In D non - ac i(t) (or /(r)) with r is identical to 
that of n, which is not necessarily given by the traditional 
KZ scaling. Does this scenario hold true in general? Be- 
low we highlight a special situation where this connection 
between the scaling of D non _ a d{t) and n, clearly breaks 
down. 

Let us probe the scaling of D n on-ad when the param- 
eter 7 of the environment ([T]) is quenched as 7 = t/r 
so that the spin chain is swept across the anisotropic 
critical point (for \h\ < 1) and the MCPs along the 
gapless Ising transition lines for \h\ = 1 (see Fig. [T|). 
Note that here one rewrites Eq. ((T|) in terms of 7 with 



J x 1 'Jy 



1, and modifies H$e to the form Hse 



-(<V 2 ) EiOf of+i ~ o'M+iVt- This represents a CSM 
in which the CS couples to the XY spin chain through 
the parameter 7. The coupling 5 therefore provides two 
channels of the temporal evolution of the environmental 
ground state with anisotropy 7 + S and 7 — 6, respec- 
tively. We recall that the problem was studied in Rcf ^22, 
from the view point of defect generation. For \h\ < 1, 
n ~ t -1 / 2 , as expected from KZ theory. For the DF, one 



numerically verified (see Fig. [2jb)). Surprise emerges for 
\h\ = 1 where one finds n ~ t -1 / 3 , a scaling that can not 
be explained in terms of traditional KZ theory Moreover, 

it was Shown that Pk = e ~^r(l+coskf / sink _ e -7TTfc 3 /2 

for k ~ 7r when h = 1. Docs this imply a scal- 
ing In D non - a d(t) ~ t 2 Jt 1 ^ for gapless quenching (see 
Fig. 1)? 

To address this question, we explore h = 1 case in 
details. Using an appropriate basis^S, one can recast the 
reduced (2 x 2) Hamiltonian Hk{t) to the form 



Hf(t) 



= 2 



(7 ± S) sin k 
h + cos k 



h + cos k 
- (7 ± 6) sin k 



(7) 



Using Eq. §4§ anc | noting that A = A8k we find that 



F k = 1 - 4(e" 7rrfc3/2 



)sin 2 (4(Sfci) 



(8) 



for the modes close to k — ir. Assuming the limit 6 — > 
and using mathematical steps identical to those employed 
in deriving Eq. (6) starting from Eq. (5), we once again 
find an exponential decay given by 



D non „ ad (t) ~ exp{-2 14 / 3 ^ 2 i 2 /(37rr)}. 



(9) 



finds that \nD. n 



1— ad 



t 2 /r 



1/2 



7 



2 r 3 / 2 , which is also 



We therefore find a clear deviation in the scaling of 
Dnon-ad (or /(t) ~^ t _1 ) from n ~ r -1 / 3 . In the present 
case, the momentum dependence of the term sin (4Skt) 
in Eq. ((5J) renders an additional r -2 / 3 factor resulting 
to a 1/r scaling of lnD„ „_ ac ;. This clearly presents a 
situation where there is no direct connection between 
n and D n on-ad- Substituting t = 77-, one finds that 

lnD non _ ad (2 14 / 3 7V(5 2 7 2 r)/37r; this is numerically 

verified as shown in Fig. [2Jb). We note that the scaling 
© can also be reproduced analytically by solving the 
Schrodinger equation with equivalent reduced Hamilto- 
nian H k in Eq. Q. 

In conclusion, we have found out the scaling of the DF 
(or /(t)) of a qubit coupled to a quantum spin chain 
which is driven across QCPs and quantum critical lines. 
We show that the scaling of the DF is given by the scaling 
of n for linear and non-linear quenching through isolated 
critical points. More importantly, our studies also reveal 
that this scenario is not universally valid. 

AD and UD acknowledge CSIR New Delhi, for finan- 
cial support. 
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